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ABSTRACT 

The Chandra X-ray Observatory obtained a 50-ks observation of the central region of M81 using the 
ACIS-S in imaging mode. The global properties of the 97 x-ray sources detected in the inner 8'. 3 x 8'. 3 
field of M81 are examined. Roughly half the sources are concentrated within the central bulge. The 
remainder are distributed throughout the disk with the brightest disk sources lying preferentially along 
spiral arms. The average hardness ratios of both bulge and disk sources are consistent with power law 
spectra of index r~1.6 indicative of a population of x-ray binaries. A group of much softer sources are 
also present. The background source-subtracted logTV-logS' distribution of the disk follows a power law 
of index —-0.5 with no change in slope over three decades in flux. The logAT-log5 distribution of the 
bulge follows a similar shape but with a steeper slope above ~ 4 x 10 37 ergs s . There is unresolved 
x-ray flux from the bulge with a radial profile similar to that of the bulge sources. This unresolved flux is 
softer than the average of the bulge sources and extrapolating the bulge log./V-log5 distribution towards 
weaker sources can only account for 20% of the unresolved flux. No strong time variability was observed 
for any source with the exception of one bright, soft source. 

Subject headings: X rays: galaxies — galaxies: individual (M81) — galaxies: luminosity function — 
galaxies: structure — galaxies: fundamental parameters 
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deferred to a future paper. 

2. OBSERVATIONS 

A 49926 second observation of a portion of the galaxy 
M81 was obtained with the Chandra Advanced CCD 
Imaging Spectrometer (ACIS, G. P. Garmire et al., in 
preparation) on 2000 May 7. The primary target, SN 
1993J, was located near the nominal aimpoint on the back- 
illuminated CCD S3 of the spectroscopy array (operating 
in imaging mode). The nucleus of M81 lies 2/79 from 
SN 1993J towards the center of S3 in this observation. 
The observation was taken in faint time exposure mode 
at 3.241 s-frame . Standard Chandra X-ray Center pro- 
cessing has applied aspect corrections and compensated 
for spacecraft dither. We selected the standard grade set 
and events in pulse invariant (PI) channels corresponding 
to ~ 0.2 to 8.0 keV for analysis. 

The majority of sources identified are within the 8'. 3 x 
8'. 3 field of S3. To avoid difficulties with the intrinsic dif- 
ferences between front- and back-illuminated devices we 
restrict our analysis to CCD S3. Due to poor cosmic ray 
rejection for the central two columns of S3, all events with 
CHIPX equal to 512 or 513 were excluded from analy- 
sis. Since these columns are dithered across the sky, their 
removal results in a slight loss of sensitivity (5%) in the 
central 16". Pixels with centers less than 1".25 from the 
readout trail of the nucleus were flagged and excluded from 
analysis as were pixels containing less than 50% sky expo- 
sure near the edge of the CCD. 

Sources were located using a minimum variance estima- 
tor that compares the data to the telescope Point Spread 
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1. INTRODUCTION 

Among the best-studied galaxies beyond the Local 
Group is the nearby Sab spiral M81 (NGC 3031). With an 
inclination of 59°, a dominant two-armed spiral pattern, 
and a well-defined central bulge, M81 is ideal for investi- 
gating formation and evolution of galaxies similar to our 
own. 

Most x-ray studies of M81 have concentrated on the 
dominant source in the field, a low-luminosity AGN at its 
core. The only extensive examination of the entire M81 
field at x-ray energies is that of Fabbiano (1988, hereafter 
F88) based on Einstein observations. Einstein detected 9 
sources in M81 including the nucleus, all with x-ray lu- 
minosities exceeding ~ 2 x 10 38 ergs s _1 , with 5 of these 
located along spiral arms. We have identified an additional 
9 sources within a ~6' radius of the nucleus from archival 
ROSAT observations. 

In this first examination of a Chandra X-ray Observa- 
tory observation of M81, we analyze the global properties 
of the 97 x-ray sources detected in the inner 8'. 3 x 8'. 3 field 

ergs- 



to a limiting 0.2-8 keV luminosity of L ~ 3 x 10 36 



s . Details of the observation and the analysis methods 
are given in § ^. Section |^ examines the spatial distribu- 
tion of sources and identifies those sources spatially coinci- 
dent with previously classified objects. The hardness ratio 
distribution is presented in § [| and the log-ZV-logS 1 distribu- 
tion in § |^. Section ^| discusses the luminosity evolution of 
sources common to the Einstein (c.1980), ROSAT (-1990) 
and present Chandra observations. The results are sum- 
marized in § ^. Detailed analysis of individual sources are 
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Function (PSF). This method is most suited to finding 
point sources that roughly match the telescope's PSF. The 
PSF was approximated by a circular gaussian with the 
same width as the measured PSF at the aimpoint and in- 
creasing quadratically with distance from the aimpoint to 
approximately match the core and off-axis broadening of 
the true PSF. We define the signal-to-noise ratio (S/N) to 
be the estimated source counts divided by the la uncer- 
tainly in this number. To set a minimum threshold, we 
considered the number of sources found per S/N interval 
down to a minimum S/N of 1.0. This distribution shows a 
clear gaussian noise peak at low S/N that indicates a sin- 
gle false detection (on S3) is to be expected for a minimum 
S/N of 3.0 and less than 0.1 false detections at a threshold 
of 3.5. We include all sources with S/N>3.0 in our anal- 
ysis with the exception of the spectral analysis (§ ||) and 
logA-logS* analysis (§ ||) for which only S/N>3.5 sources 
are considered. 

A second pass, designed to mitigate source confusion, 
was made through the data. For this purpose, all pixels 
within 2a of each source position were flagged and the 
source-finding process repeated. In this way, the variance 
in the local background due to nearby sources is greatly 
reduced and otherwise confused sources isolated. An ad- 
ditional 6 sources, 4 with S/N>3.5, were detected. 

The final source list contained 81 sources with S/N>3.5 
with an additional 16 with 3.0<S/N<3.5. These source 
positions are identified in the ACIS data shown in fig- 
ure la and superimposed on the second generation digital 
sky survery (DSS) image in figure lb. 

Small spatial regions centered on each source were then 
selected for further analysis. The size of these spatial re- 
gions also varies quadratically with off-axis position and 
encompasses roughly the 95% encircled energy radii of the 
PSF at 1.5 keV. No correction for the low background, 
~0.04 cts-pixel -1 , is made. For spectral analysis (§ Q), the 
total counts within each source region were binned into soft 
(0.2 - 2.0 keV) and hard (2.0 - 8.0 keV) bands. Higher res- 
olution spectra, at 16 PI channel binning out to ~10 keV, 
were examined for spurious effects such as anomolously 
high background. No anomolies were encountered. 

Finally, we distinguish between the 'bulge' and 'disk' 
defined here as those regions internal and external, respec- 
tively, to a 4/70 x 2/35 ellipse centered on the nucleus with 
major axis position angle 149° (cf. figure 1). This ellipse 
corresponds to a circle in the plane of the galaxy of radius 
2/35 (^2.5 kpc). There is extended, unresolved x-ray flux 
in the bulge that is above background and is not due to the 
PSF wings of the nucleus nor to other point sources. This 
component will be referred to as 'excess bulge emission.' 

3. MORPHOLOGY 

The bulge contains 41 sources at a surface density of 5.4 
sources-arcmin -2 and the disk contains 56 sources (on S3) 
at a density of ~1.0 source-arcmin -2 . Among the latter, 
21 are located within ±400 pc of spiral arms (as defined 
by Matonick & Fesen (1997) by tracing the Ha emission 
peak, their figure 31b). The surface density of spiral arm 
sources in our sample is slightly less than, but statisti- 
cally consistent with, the surface density of disk sources 
not associated with spiral arms. 

The surface density distribution of bulge sources, 
inclination-corrected and centered on the nucleus, can be 



described by an exponential radial distribution, cr(r) oc 
e~ r / r °, with r Q ~ 0.92 ± 0.26 kpc. After subtracting 
the background component and the contribution from the 
wings of the PSF of the nucleus, the surface brightness 
profile of the excess bulge emission can also be described 
by an exponential with scale length r Q ~ 0.83±0.20 kpc for 
0/8 < r < 5/5. The excess bulge emission is much steeper 
at smaller radii though a residual contribution from the 
wings of the nuclear PSF cannot be definitively excluded. 
Beyond r ~ 1 kpc, these profiles are roughly consistent 
with the V (Georgiev & Getov 1991), g (Frei et al. 1996), 
and B and I (Elmegreen & Elmegreen 1984) profiles sug- 
gesting the bulge x-ray emission follows the distribution of 
starlight. 

A search of the literature for known objects spatially 
coincident with Chandra sources found no matches within 
the bulge, with the exception of the nucleus. Among the 
disk sources, 5 are located within 3" of supernova rem- 
nants (SNR; Matonick & Fesen 1997) including SN 1993J 
and the bright Einstein source X-6 (F88). All 5 SNR can- 
didates lie on spiral arms as expected for young, massive 
star, Type II supernova progenitors. Source X-6 is a point 
source located within a 90-pc-diameter SNR (Matonick & 
Fesen 1997). The high luminosity and hard spectrum (§ [|) 
of X-6 is uncharacteristic of SNRs. 

Six x-ray sources, all located within spiral arms, are 
within 10" of known HII regions, half of which are giant 
radio HII regions (Kaufman et al. 1987). One of these x- 
ray sources, though weak, has the hardest spectrum of the 
entire source sample (§§]). 

No x-ray sources were found to coincide with known 
globular clusters, novae, or stars. Roughly 10% of globu- 
lar clusters are expected to contain accreting neutron stars 
detectable above our threshold and galaxies like M81 typ- 
ically contain of order 200 globular clusters (Perelmuter & 
Racine 1995). However, there are only 25 confirmed glob- 
ular clusters in M81 (Perelmuter, Brodie & Huchra, 1995) 
and only 4 within the S3 field. 

The nova rate in M81 is probably of order 20 year -1 
(A. W. Shaffer, private communication). The x-ray sig- 
natures and light curves of novae are only sparsely doc- 
umented (Ogelman, Krautte & Beuermann 1987) making 
novae difficult to identify in our source sample. 

There are two sources, Einstein source X-2 and a 
ROSAT-detected object, which are near stars or star 
clusters. X-2 is l."8 distant from a star cluster with 
a spectrum consistent with G-type stars (Sholukhova et 
al. 1998). Zickgraf & Humphreys (1991) list the opti- 
cal object as a probable foreground F-type dwarf. The 
ROSAT source is coincident with an undocumented ~18- 
magnitude object visible in the Digital Sky Survey image. 
Interestingly, both x-ray sources are transient and have 
been much brighter in the past (§ [|). 

No known background or foreground objects were found 
coincident with the S3 Chandra sources. Up to ^25 
background objects are expected in the S3 field above 
S/N = 3.5 based upon our analysis of the calibration field 
CRSS J0030.5+2618 (see also Brandt et al. 2000). This 
estimate ignores any obscurtion by M81. Based on results 
of the Einstein stellar survey (Topka et al. 1982), we es- 
timate roughly 0.4 F- and G-stars could be detectable in 
our field but these should all be extremely soft (Maggio 
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et al. 1987, Hodgkin & Pye 1994, Schmitt 1997). A soft 
spectrum was observed from 3 catalogued stars detected 
onCCDS2. 

4. HARDNESS RATIO 

The total source counts in the hard x-ray band (2.0 - 8.0 
keV) is shown in figure 2 against the soft (0.2 - 2.0 keV) 
band for the subset of sources with S/N>3.5 excluding the 
nucleus. There is no adjustment for pileup which tends to 
harden spectra of the brightest sources. 

Also shown are the hardness ratios of the sum of all the 
bulge sources, and of all the disk sources, the hardness ra- 
tio of the excess bulge emission (background and nuclear 
PSF subtracted), and of the typical background (scaled to 
a 2' x 2' area). The nucleus and the bright soft source are 
omitted from the summed bulge point, and X-6 is omitted 
from the summed disk point as these sources are clearly 
spectroscopically atypical. 

The line shown corresponds to an absorbed power law 
of spectral index T = 1.6, typical of x-ray binaries, with 
the Galactic absorbing column density of Njj = 4x 10 20 
cm~ 2 (Stark et al. 1992). Though indicative and useful to 
guide the eye, this canonical spectral shape is not a model 
fit result. Most sources fall near this curve. Notable ex- 
ceptions are the soft sources at the lower right of figure 2. 
The softest bulge sources have no known associations but 
the 3 softest disk sources all lie on spiral arms. The softest 
source is also time- variable (§ ||) and is the third brightest 
source in the entire sample. The hardness ratios for the 
summed bulge and summed disk spectra are also typical 
of x-ray binaries. The excess bulge emission is softer in- 
dicating some of the excess emission is from truly diffuse, 
hot, interstellar gas and not entirely from unresolved point 
sources. 

5. LOG7V-LOGS 

The number of sources above a limiting brightness, 
N(>S), is shown against the 0.2 - 8.0 keV count rate, 
S, in figure 3. The nucleus is excluded as are sources with 
S/N<3.5. Background sources were taken from the cali- 
bration field CRSS J0030.5+2618 (§ |), by applying our 
same analysis methods, and have been subtracted after 
scaling to the fractional areas of the bulge and disk. 

The curve shown in figure 3 represents the best power 
law fit to the disk distribution, N = OASS' 050 , over 
the entire range of S. Similarly, N — 0.44S' -0 ' 57 for the 
bulge source distribution for S<0.004 count-s -1 . There is 
a break in the slope of the distribution of bulge sources 
above this point but no such break in the disk source dis- 
tribution. Of the 10 disk sources with 5>0.004 count-s -1 , 
7 are coincident with spiral arms. 

The luminosity of each source can be estimated by as- 
suming the r = 1.6 power law spectral model of § H and 
a distance of 3.6 Mpc to M81 (Freedman et al. 1994). 
For our 50-ks observation, therefore, S = 0.004 count- 
s' 1 in the 0.2 - 8.0 keV band corresponds to a luminosity 
Lx = 3.7 x 10 37 ergs s _1 and the faintest source in the 
field corresponds to a luminosity of ~ 3 x 10 36 ergs s _1 . 

Excluding the nucleus, the total bulge luminosity is 
Lx ~ 2.4 x 10 39 erg s _1 of which 36% is excess (unre- 
solved) emission. The total luminosity of the disk sources 
is Lx ~ 3.9 x 10 39 erg s _1 . The luminosity of the nucleus 



is Lx ~ 4 x 10 40 erg s _1 based on spectral fits to the 
trailed image. This is within the ASCA-observed range of 
luminosities (Ishisaki et al. 1996) and comparable to the 
BcppoSAX observed luminosity (Pellegrini et al. 2000). 
The nucleus contributes approximately 86% of the total 
luminosity in the 8'. 3 x 8'. 3 S3 field of view. 

6. TIME EVOLUTION 

A simple test for time variability was made by binning 
the light curves of all sources on 1000, 2000, and 4000 
second intervals. Applying a \ 2 test 101 the hypothesis 
of consistency with a mean value found only one clearly 
significant deviation on all three timescales. This source 
is also the softest source on S3 and the third brightest. 
It is present in at least 3 of 6 ROSAT HRI observations 
spanning 1993 - 1998 but is too close to the nucleus to be 
identified in any other previous x-ray observation. 

M81 has been observed in x-rays at moderate spatial res- 
olution by Einstein in 1979 and by ROSAT over the period 
1991-1998. There are 6 Einstein sources (F88) within the 
S3 field of view and another 4 ROSAT-detected sources 
according to our analysis. Chandra resolves six ROSAT 
source regions into two or more sources including Einstein 
sources X-7 and X-10, and, of course, the nucleus. 

Long-term variability has been detected in the Chandra 
variable source, the nucleus, X-2 and the ROSAT-detected 
source coincident with an undocumented star-like object 
(§ ||). Both of the latter sources are moderately weak in 
the present observations (~3 and ~ 6 x 10 37 ergs-s -1 , re- 
spectively) but have been much brighter in the past. 

The location of another Einstein source, X-12, places it 
on the eastern edge of S3. There are no Chandra sources 
within the portion of the 45" error circle of this IPC source 
that falls on S3 and it is not detectable in the ROSAT ob- 
servations. 

7. SUMMARY 

The global properties of the Chandra sources identified 
in the central 8'. 3 x 8'. 3 region of M81 have been exam- 
ined. There is a high density of sources located in the bulge 
of the galaxy. These sources, and the excess bulge x-ray 
emission, follow an exponential radial distribution about 
the nucleus with an ~1 kpc scale length. This is similar 
to that of optical light and suggests the excess bulge emis- 
sion is from unresolved point sources and that the point 
sources trace the old stellar bulge population. However, 
extrapolating the bulge logTV-log.!? distribution towards 
weaker sources predicts ^0.022 source counts-s from 
sources <0.001 counts-s _1 -source _1 . Thus, weak sources 
can only account for 20% of the total excess bulge emis- 
sion (^0.092 counts-s -1 ). There remains a significant flux 
in the bulge unaccounted for by unresolved sources, par- 
ticularly at lower energies, indicative of a truly diffuse, 
hot interstellar gas. If all the excess bulge emission arises 
from a diffuse thermal bremsstrahlung at kT ~ 0.2 keV, 
uniformly distributed throughout the bulge, then there is 
- 4 x 10 6 M of hot gas. 

The average hardness ratio of the Chandra sources can 
be described by a power law of index T = 1.6 with a Galac- 
tic absorbing column density of Njj — 4 x 10 20 cm . 
This is typical of x-ray binaries. Assuming thie canoni- 
cal spectrum applies to all sources, the background source 
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subtracted logTV-logS distribution of the disk sources fol- 
lows an N oc profile extending to the most lumi- 
nous sources at Lx > 10 39 ergs-s -1 . The distribution of 
bulge sources, however, shows a break at L x ~4x 10 37 
ergs-s _1 . Thus the bulge sources are predominantly low- 
mass x-ray binaries, consistent with an old bulge stellar 
population, while the disk has an additional population of 
much brighter sources. Seven of the 10 brightest Chandra 
disk sources lie within spiral arms, consistent with previ- 
ous Einstein observations (F88), but the entire sample of 
disk sources, down to the limiting Chandra luminosity of 
~ 3 x 10 36 ergs-s _1 , are distributed throughout the disk 



without preference to spiral arms. Some of these brighter 
sources may be black hole candidates with high-mass com- 
panions. 

We thank Martin Weisskopf for deriving the expression 
for the uncertainty in the source flux used in the source- 
finding algorithm. The digital sky survey image used in 
Figure lb is from the "Palomar Observatory - Space Tele- 
scope Science Institute Digital Sky Survey" of the northern 
sky, based on scans of the Second Palomar Sky Survey, and 
was produced under NASA Contract NAS5-2555. 
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FIGURE CAPTION 

Fig. 1. — Left: Chandra ACIS-S image of the inner region of M81 in J2000 coordinates. North is up and east is to the right. Each sky pixel 
is 0.492". The aimpoint is near SN 1993J. The readout trail, caused by photons from the bright nucleus hitting the detector during readout, 
has been removed. Chandra sources are marked as light blue (S/N>3.5) or magenta (3<S/N<3.5) circles. Einstein sources (Fabbiano 1988) 
are marked by red circles encompassing the approximate Einstein positional uncertainties. The 4'. 7 X 2'. 35 ellipse denotes the boundary 
between 'bulge' and 'disk' as used in the text. Right: 73-band digital sky survey image of M81 with the locations of sources identified as in 
left panel. 

Fig. 2. — Hardness ratio of the 80 sources with S/N>3.5, excluding the nucleus, with v^V errors. Also included are the sum of all 
bulge sources (B) excluding the bright soft source and the nucleus; the sum of all disk sources (D) excluding Einstein source X-6; the excess 
(unresolved) bulge emission (o); and the average instrument background (*) over a 2' X 2' region. The line denotes the hardness ratio of a 
typical x-ray binary; a Galactically-absorbed power law of index T = 1.6. 

Fig. 3. — Background-source-subtracted logAf-logS 1 for bulge and disk sources. The background log/V-logS distribution was taken from 
Chandra calibration observation of CRSS J0030. 5+2618 (sec also Brandt et al. 2000) using the same analysis procedures used in this study 
of M81. The line represents the power law fit to the disk distribution, TV = O^SS 1-0 ' 50 . 
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